Two-pion correlations in Au+Au collisions at 10.8 GeV/c per nucleon 
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Two-particle correlation functions for positive and neg- 
ative pions have been measured in Au-|-Au collisions at 
10.8 GeV/c per nucleon. The data were analyzed using one- 
and three-dimensional correlation functions. From the results 
of the three-dimensional fit the phase space density of pions 
was calculated. It is consistent with local thermal equilib- 
rium. 

Bose-Einstein correlations of identical pions can be 
used to obtain information about the space-tinie- 
niomentum distribution S{r,t,p) of pions at freeze-out 
(pion source) in a nuclear collision (for a review see 
Q). The two-pion correlation function C(pi,p2), de- 
fined as the ratio between the two-particle density and 
the product of single particle densities, shows a peak at 
q = P2 — Pi =0. The width of this peak is inversely 
proportional to the size of the pion source. Analysis 
of C as a function of different components of the rela- 
tive momentum yields information about the source di- 
mensions in different directions. The relation between 
S and C is, under realistic conditions, complicated by 
effects like three- (and more) particle correlations [gj-Q, 
long-lived resonance decays [^|-0, and the distortion of 
the single-particle spectra by the two-particle correla- 
tions H,^. The correlation must also be corrected for 
the Coulomb interaction between the two pions ]l0|-fi2[ . 

In this Letter we present results of the correlation 
analysis of pions produced in central Au-t-Au collisions 
at 10.8 GeV/c per nucleon. The data were taken in 
Fall 1993 at the AGS. The central trigger at the level 
of 10% of the geometrical cross section was used. The 
apparatus allowed a simultaneous measurement of posi- 
tive and negative pions with a momentum resolution of 
Ap/p « 3%. With two field polarities used, the overall 
acceptances for positive and negative pion pairs were sim- 
ilar: 2 < y < 4 and < pt < 0.5 GeV/c with (y) = 3.1 



and (pt) — 0.1 GeV/c. The analyzed data sample con- 
sists of one million central events and the total number of 
analyzed 7r+7r+, tt~tt~ , and tt^tt" pairs is 130 k, 210 k, 
and 340 k, respectively. The particle identification (PID) 
quality was tested by varying the PID cuts in order to 
deliberately accept background particles. No significant 
influence on the results was observed implying that par- 
ticle misidentification does not contribute significantly to 
the overall systematic uncertainty. A more detailed de- 
scription of the experiment and of the data analysis can 
be found in |13| . 

Experimentally, the correlation function C(q) is de- 
fined as the number of pion pairs in a q-bin (signal) di- 
vided by the number of such pairs obtained by event 
mixing (background). The variables ^outiQside, and (7iong, 
used in the three-dimensional 'out-side-long' analysis, 
are the components of q in the beam rapidity frame 
2/anai = ybcam = 3.14. Here giong is the component par- 
allel to the beam, qside is perpendicular to the beam and 
to the average pair momentum, and gout is perpendic- 
ular to giong and gsidc Mjfsl- In our one-dimensional 
analysis we analyzed two-pion correlations as a func- 
tion of q — |q| calculated in the pair cm. frame. For 
equal mass particles q is equal to Qinv which is defined 
as V(P2-Pi)'-(^2-£^i)2. 

The background pairs were weighted with the Coulomb 
correction factor to remove the effect of the Coulomb 
interaction between the two pions from the correlation. It 
is calculated by taking the square of the non-relativistic 
wave function describing a particle in a Coulomb field 
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where F is the confluent hypergeometric function, k is 



the asymptotic momentum of a pion in the pair c.m.s. 
{k = Qinv/2), and r is the relative distance of the two 
pions at freeze-out. The relative velocity enters through 
77 = ZiZ2m.^ca/ Qi-a;^ where ttt-tt is the pion mass and a 
is the fine structure constant. For pions emitted from a 
point-like source F is equal to unity and the Coulomb 
correction becomes the well known Gamow factor [ p^|Jl6| 
(dot-dashed line in Fig. 1). For a finite size source, the 
Coulomb correction factor can be obtained numerically 
by averaging i?(k, r) over r. The correction H calculated 
this way for a Gaussian source with ax = Uy = Gz = 
5 fm and at — Q is shown as a dashed curve in Fig. 1. 
In addition, since it would be incorrect to divide data 
which have been measured with a finite resolution Ap 
by a correction calculated assuming Ap=0, we folded H 
with the momentum resolution of the spectrometer. The 
resulting correction H\p is shown as a solid line in Fig. 1. 
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FIG. 1. Coulomb correction factor for two-pion correla- 
tion functions. The dot-dashed line G represents a correc- 
tion assuming a point-like pion source (Gamow factor). The 
dashed line H was calculated for a Gaussian source with a 
realistic size. The same curve folded with the spectrometer 
momentum resolution is shown as a solid line. This is the 
correction used in the data analysis. 

Since the tt+tt" correlation is dominated by the mutual 
Coulomb interaction, we can use it to test the quality of 
the Coulomb correction. This correlation function, un- 
corrected, corrected by the Gamow, and by H^p is shown 
in the right hand panels of Fig. 2. The properly corrected 
correlation should be equal to unity. The Gamow correc- 
tion is obviously inappropriate. This is in agreement with 
several theoretical studies |ll|,|l8| , |l9| as well as previous 
experimental observations po| , pi| . The by i?Ap corrected 
7r+7r~ correlation is close to unity. We therefore use the 
same method to correct the tt+tt^ and 7r~7r^ correlation 
functions. 
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FIG. 2. One-dimensional tx^-k^ ,-k~-k~ , and tt^tt" corre- 
lation functions. Presented are from top to bottom: raw and 
Gamow corrected correlations, and correlations corrected by 
H/^p. The solid line is the Gaussian fit. 

The Gamow correction would lead to a significantly 
different shape of the corrected correlation (see Fig. 2 
and Table I). 

In order to test possible distortions of the correlation 
function introduced by the apparatus and the data anal- 
ysis we performed a Monte Carlo simulation in which 
we generated events with realistic particle multiplicities 
and momenta, and propagated the particles through our 
setup using GEANT |g^. The simulated detector re- 
sponse was used as an input to the data analysis code. 
The obtained correlation functions were equal to unity 
as expected for a simulation without two-particle inter- 
actions. Subsequently, we generated events in which posi- 
tive pions were correlated. The obtained correlation func- 
tions were fitted and the fit results were compared to the 
parameters of the original correlation. The observed 3- 
10% and 20% reductions of the radius parameters and the 
correlation strength (see below for definitions), respec- 
tively, are consistent with the broadening of the correla- 
tion peak expected from the finite momentum resolution 
of the spectrometer. 

The experimental. Coulomb corrected one-dimensional 
7r^7r+ and 7r~7r~ correlations were parametrized by 
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with h = 1. The results of a maximum hkehhood fit 
are shown as the sohd line in Fig. 2. The fit parameters 
are given in Table I. Fitting Coulomb uncorrected cor- 
relations gives the same R but a lower A value of 0.32. 
We estimate the uncertainty of A because of the assump- 
tions entering the Coulomb correction to be half of the 
difference, i.e. 18%. 

In the 'out-side-long' analysis, which was performed 
in the beam rapidity frame j/anai — 3.14, the three- 
dimensional correlation functions were fitted by 



C(9out, 9sido, fflong) 



1 + Aexp(— i?out Sout — ^sido (Zsidc 
— ^long fflong — 2|i?ol|-Rol'?out<Zlong) • 



The projections of the correlation and of the fit, obtained 
by averaging over the two other components of q in the 
range from -50 MeV/c to 50 MeV/c with a weighting fac- 
tor equal to the number of background pairs, are shown 
in Fig. 3. 
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Figure 3 
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FIG. 3. Projections of the three-dimensional tt'^tv^ and 
TT~n~ correlation functions. Solid lines represent projections 
of the fit. 

The particles were labeled such that giong>0. The 
asymmetric shape of the projection C(qout), well repro- 
duced by the fit, is caused by the acceptance and the 
finite width of the projection window. The parameters 
obtained from the fit are given in Table IL The values 
were corrected for the momentum resolution of the spec- 
trometer using the results of the Monte Carlo simulation. 

Similarly as in the one-dimensional fit, the tt^ and tt" 
sources apparently differ in size. We attribute this to 
the distortion of the correlation function caused by the 
Coulomb interaction between each of the pious and the 



rest of the system [Q . 

The parameter Rgidc is most directly related to the 
transverse size of the pion source. In order to compare it 
to the projectile size we first translated the r.m.s. radius 
of a gold nucleus i?rms=5.3±G.l fm Eq] to an equivalent 
Gaussian radius Rg — l/-\/3i?rms = 3.08 ± 0.06 fm. The 
ratio Rsidc/ Rg is 1.3±0.3 for tt"*" and 1.8±0.2 for 7r~. 
This increase of the transverse size, combined with the 
recent finding of transverse flow velocities from analysis 
of particle spectra p4,M, may be used to restrict the 
expansion dynamics. 
+ It has been shown recently that the two-meson correla- 
(35ion function is related to the meson phase space density 
(defined as the number of particles per unit volume h^ 
in six-dimensional phase space) at freeze-out |26 . Using 
the results from |2q] we obtained the following formula: 
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where (/)p is the phase space density of pious 
with momentum p, averaged over the spatial coordi- 
nates, (fin/(fip is the differential pion multiplicity, and 
A, i?out, -Rside, -Riong, and Rqi are the numbers obtained 
from the 'out-side-long' fit to the correlation function. 
This phase space density can be compared to a theo- 
retical prediction assuming local equilibrium, i.e. to the 
Bose-Einstein distribution function. However, pious from 
long-lived resonances, which result in A < 1 [Q,^, de- 
crease the mean phase space density. This component of 
the source function is not thermalized and should be left 
out from the comparison. Since for a case without long- 
lived resonances A = 1, and since the actual fraction of 
pious coining from the core is vA, the following formula 
describes the average phase space density in the core: 
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Taking the correlation parameters from Table II (with 
the error matrix properly taken into account) and us- 
ing experimental cPn/cPp from Ref. |g^, we calculated 
(Z)™"" for several points in the {pt,y) plane (Table III). 
We assume that the correlation parameters do not change 
significantly between the points and thus the entire mo- 
mentum dependence comes from (Pn/dPp; a more correct 
approach would require a separate correlation analysis for 
every {j>t,y) bin. 

In the absence of tranverse flow and neglecting spectral 
changes due to resonance decay, the phase space density 
of pious in thermal equilibrium would follow the Bose- 
Einstcin distribution (for each pion species): 

^ ^P^"exp(miCosh(y-2/s)/T)-l ^^^ 

with the pion transverse mass rrit and rapidity y, and 
the source rapidity ^5. Since transverse momentum spec- 
tra of pious and protons at these rapidities have similar 



slopes 1^,^ , the neglect of transverse flow seems to be 
appropriate. Nevertheless, to estimate an upper limit for 
the influence of flow effects we found by Monte Carlo 
simulation that transverse flow with /3=0.3 results in a 
phase space density lower by 20% than the one given by 
Eq.(p|). Deviations of the spectral shape from thermal 
distribution due to A decays could be handled best if the 
analysis was done in bins of pt and y. This is planned 
for the high-statistics data sets taken in the '94 and '95 
runs of E877. 

Under these assumptions, we can test the thermaliza- 
tion of pions by comparing the experimental phase space 
density to that given by Eq. (|6|) . The comparison is sim- 
plifled by the fact that the experimental differential mul- 
tiplicity can also be parametrized by the Bose-Einstein 
function: 



d^n _ A 

(Pp exp(mt/Tcff) - 1 
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with two fit parameters: effective transverse tempera- 
ture Tcff and normalization A. We assume that Tcff = 
T/ cosh(2/ — ys)- Consequently, the comparison between 
(/)p°'^° and /^^(p) is reduced to the comparison between 
the parameters of the two-pion correlation function and 
the normalization of the pion spectra A: 
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Since the p^-spectra have a somewhat concave shape, A 
depends on the range of the fit: low (high) pt's yield high 
(low) A. Using A obtained from a fit to the entire mea- 
sured transverse momentum spectra < pt < 0.6 GeV/c 
in the rapidity range 3.0-3.3 and the correlation param- 
eters from Table II, we evaluated the experimental-to- 
thermal density ratio of 0.97±0.23 and 1.02±0.18 for 7r+ 
and 7r~, respectively. Thus the experimental pion phase 
space density is consistent with the presence of local equi- 
librium. 

Summarizing, the multi-dimensional pion source pa- 
rameters have been determined in the Au-|-Au system at 
the AGS. The observed source is significantly larger than 
a gold nucleus. The extracted pion phase space density 
at freeze-out indicates local equilibrium. 
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TABLE I. Parameters of the fit to one-dimensional cor- 
relation functions C(Qinv) = 1 + Aexp(— i? ^Qinv^). For the 
purpose of comparison we also quote the results obtained with 
the Gamow correction (denoted as G). 



with H^p correction 
A R (fm) 



with G correction 
A R (fm) 



0.49±0.04 
0.51±0.03 



5.5±0.3 
6.2±0.2 



0.56±0.04 
0.62±0.03 



5.1±0.2 
5.9±0.2 



" Present address: Gesellschaft fur Schwerionenforschung, 

Darmstadt, Germany. 
'' Present address: University of Heidelberg, Germany. 



TABLE II. Parameters of the fit to three-dimensional cor- 
relation functions C(gout,<7sidc, giong) = 



1 + Aexp(-_Rout gout 



jtside Qsidc -r Jtlong <?long 



— 2|_Roi|fioigout9iong)- The fit results are corrected for mo- 
mentum resolution. That is why A's difi'er from the values in 
Table I. 



A 



+7r+ 0.62±0.06 
-yr" 0.62±0.05 



i?out(fm) i?side(fm) i?iong(fm) i?oi(fm) 



5.8±0.5 
6.5±0.5 



3.9±0.8 
5.6±0.7 



5.5±0.4 
5.8±0.4 



2.4±0.7 
3.7±0.8 



TABLE III. Pion phase space density at freeze-out as a 
function of pt and y. The statistical errors from cPn/d^p 
are about 2%. The systematic error from (Pn/d^p, combined 
with the error from the correlation analysis, is 24% for 7r+ 
and 18% for 7r~. (This error acts on all numbers in the Table 
collectively.) 
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